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INTRODUCTION 


Voids,  cracks,  or  fissures  in  soLid  propellant  grains  have  been 
responsible  for  the  failure  of  many  rocket  motors  to  meet  acceptance 
specifications.^  These  defects  may  originate  as  a result  of  the  manu- 
facturing process,  damage  received  during  handling  and  transportation, 
and  chemical  aging  or  thermal  stresses  incurred  during  storage.  However, 
even  though  the  propellant  grain  is  crack-free  before  firing,  there  is  a 
possibility  that  cracks  can  form  in  regions  of  large  stress  concentration 
during  rocket  motor  operation.  Combustion  within  a crack  can  lead  to  a 
pressure  build-up  sufficient  to  cause  further  propagation  of  the  crack, 
resulting  in  an  unexpected  Increase  in  the  surface  available  for  burning. 
Thus,  there  has  been  much  concern  among  rocket  motor  designers  as  to  the 
extent  that  crack  combustion  may  contribute  to  catastrophic  failure. 

The  ignition  and  combustion  processes  in  solid  propellant  cracks 
generally  involve: 

1.  Penetration  of  rocket  chamber  combustion  gases  into  the  crack 
cavity . . 

2.  Convective  heating  of  the  propellant  crack  surface  to  ignition 
and  transient  conduction  within  the  solid  propellant  . 

3.  Pressure  wave  phenomena  in  the  longitudinal  direction  of  the 
crack  • 

4.  Flame  spreading  along  the  crack. 

5.  Cavity  pressurization  due  to  gasification  of  the  solid  propel- 
lant . 


^Raytheon  Company.  Stopoge  Reliability  of  Missile  Material  Progpom, 
Solid  Ppopellant  Rocket  Motor  Analysis^  by  D.  F.  Malik.  Huntsville, 

Ala.,  May  1976.  (Report  No.  LC-76-OR1  (also  DDC  Report  AD-A026105), 
publication  UNCLASSIFIED.) 
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6.  Flow  reversal  caused  by  cavity  pressures  in  excess  of  that  in  | 

the  rocket  chamber  . ' 

! 

7.  Deformation  of  the  crack  geometry  due  to  burning  and  pressure  I 

loading  . 1 

8.  Propagation  of  the  crack  tip  due  to  regression  and  mechanical  1 

fracture.  i 

] 

In  addition  to  these  interrelated  phenomena,  there  are  many  other  impor-  ] 

tant  factors  which  can  influence  the  crack  combustion  process.  These  | 

factors  Include: 

i 

1.  Initial  geometry  of  the  crack 

i 

2.  Variation  of  rocket  chamber  conditions  i 

j 

3.  Erosive  burning  effects  due  to  high  gas  velocities  '• 

\ 

4.  Dynamic  burning  effects  due  to  rapid  pressure  transients 

5.  Physlochemlcal  properties  of  the  propellant 

6.  Initial  temperature  of  the  propellant. 

Although  Imperfections  in  solid  propellant  grains  may  be  distrl-  I 

buted  in  an  irregular  Intricate  network,  a rational  starting-point  for  | 

understanding  the  complicated  interaction  of  the  above-mentioned  phe-  ] 

I 

nomena  is  to  consider  a simple  ordered  system — a single  Isolated  crack.  } 

For  such  a system,  it  is  desirable  to  establish  a criterion  for  stable  ' 

combustion  without  catastrophic  crack  propagation;  the  delineation  of  a 
safe  domain  of  governing  parameters  is  the  ultimate  goal.  But  before  it  j 

can  be  achieved,  a realistic  theoretical  model  describing  the  transient  ] 

combustion  processes  inside  a propellant  crack  has  to  be  established  and  j 

experimentally  verified.  | 

The  current  research  program  was  undertaken  to  experimentally  and  i 

theoretically  investigate  the  complex  problem  of  transient  combustion  J 

within  a propellant  crack.  The  objectives  of  the  work  for  this  report- 
ing period  are: 

1.  To  establish  a theoretical  model  for  predicting  the  rate  of 
flame  propagation  inside  a single  isolated  crack  of  variable 
geometry. 
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2.  To  predict  the  pressure  distribution  and  pressurization  rate 
inside  a crack,  leading  to  a basis  for  calculating  overall 
stress  loading  in  the  crack. 

3.  To  parametrically  investigate  the  Influence  of  crack  geometry, 
main  chamber  pressurization  rate,  and  other  physlochemlcal 
parameters  on  the  gasification  rate  inside  the  propellant 
cavity. 

4.  To  design  and  test  an  experimental  combustion  chamber  for 
verification  of  the  theoretical  model  and  for  establishment  of 
a fundamental  data  base  in  crack  combustion. 

This  report  has  been  written  to  document  the  progress  gained  toward 
achieving  a better  understanding  of  transient  crack  combustion.  The 
Background  section  reviews  some  of  the  previous  analytical  and  experi- 
mental studies  in  crack  combustion,  and  "Theoretical  Modeling"  describes 
the  boundary  conditions  and  numerical  solution  procedures.  In  "Experi- 
mental Apparatus  and  Procedure,"  the  experimental  combustion  chamber  and 
data  acquisition  system  are  discussed.  "Discussion  of  Results"  describes 
the  theoretical  results  of  a parameter  study  along  with  a comparison  of 
theory  with  limited  experimental  data.  Finally,  the  progress  achieved 
and  the  conclusions  drawn  from  the  results  obtained  are  summarized  in 
the  "Summary  of  Progress  and  Conclusions." 


BACKGROUND 


Although  the  overall  process  of  transient  crack  combustion  is 
extremely  complicated,  various  Investigations  have  been  conducted  to 
study  certain  particular  processes  described  earlier. 

Taylor  conducted  experimental  tests  to  study  the  convective  burning 
of  porous  propellants  with  closed-  and  open-end  boundary  conditions.^ 

He  observed  a critical  pressure  above  which  the  hot  gas  can  penetrate 
into  the  porous  propellant  and  significantly  Increase  the  regression 
rate  of  the  charge.  His  experiments  revealed  that  the  flame  front 
decelerated  in  a region  near  the  closed  end  of  the  propellant  charge. 


^J.  W.  Taylor.  "The  Burning  of  Secondary  Explosive  Powders  by  a 
Convective  Mechanism,"  Trans  Farad,  Soo,  Vol.  58,  1962,  p.  561. 
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Prentice  studied  the  flame  spreading  phenomena  In  transient  solid 
propellant  cracks  using  high  speed  motion  pictures  and  fine  thermo- 
couples.^ Difficulties  were  encountered  In  using  fine  thermocouples  due 
to  the  long  response  time  of  the  thermocouple,  uncertainty  In  the  pre- 
cise location  of  the  thermocouple,  and  also  due  to  the  gaseous  explo- 
sion. For  closed-end  cracks,  with  small  diameters,  he  observed  that  the 
flame  cannot  propagate  Into  the  crack  when  the  bomb  pressure  Is  low. 

Bobolev,  et  al. , studied  the  mechanism  by  which  combustion  products 
penetrate  Into  the  propellant  crack.**  Two  mechanisms  were  proposed: 

(1)  forced  penetration  of  gas  due  to  a higher  external  pressure,  fol- 
lowed by  (2)  "spontaneous  penetration"  of  product  gases  from  the  crack 
surface  Into  the  unburned  region  of  the  propellant  crack. 

Belyaev,  et  al,  showed  that  burning  of  propellant  Inside  a narrow 
pore  may  lead  to  an  excess  pressure  build-up.^  Their  experimental 
results  have  Indicated  that  the  erosive  burning  effect  cannot  be  neg- 
lected In  the  crack  combustion  study. 

In  a later  study,  Belyaev,  et  al.,  made  a series  of  experimental 
tests  to  study  the  dependence  of  flame  spreading  rate  on  crack  geometry, 
propellant  properties,  boundary  conditions,  and  combustion  chamber 
pressures.^  Part  of  their  experiments  were  conducted  In  the  atmosphere 
and  In  a "constant  pressure"  Crawford  vessel  where  the  Initial  pressure 
In  the  crack  was  equal  to  the  bulk  pressure  In  the  vessel.  They  found 
that  both  the  flame-spreading  rate  and  pressure  gradient  along  a crack 
are  strong  functions  of  crack  width.  For  a closed-end  crack,  they 
observed  that  the  flame-spreading  rate  Increases  Initially  until  reach- 
ing a constant  speed  and  then  decelerates  near  the  tip  of  the  crack. 


^Naval  Ordnance  Test  Station.  Flashdown  in  Solid  Propellants,  by 
J.  L.  Prentice.  China  Lake,  Calif.,  NOTS,  December  1962.  (NOTS  TP 
3009,  publication  UNCLASSIFIED.) 

K.  Bobolev,  A.  D.  Margolin,  and  S.  V.  Chulko.  "The  Mechanism 
by  Which  Combustion  Products  Penetrate  Into  the  Pores  of  a Charge  of 
Explosive  Material,"  Doklady  Akademii  Nauk,  USSR,  Vol.  162,  No.  2,  May 
1965,  pp.  388-391. 

^A.  F.  Belyaev,  et  al.  "Development  of  Combustion  In  an  Isolated 
Pore,"  Combustion,  Explosion  & Shook  Waves,  Vol.  5,  No.  1,  Jan-Mar  1969, 
pp.  4-9. 

®A.  F.  Belyaev,  et  al.  "Development  of  Burning  In  a Single  Pore," 
Transition  of  Combustion  of  Condensed  Systems  to  Detonation,  Chapter  5, 
Part  A,  Section  22,  Science  Publisher,  1973,  pp.  115-134. 
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Cherepanov  stated  that  In  a sufficiently  narrow  and  long  cavity,  as 
a result  of  the  Impeded  gas  flow,  the  pressure  reaches  such  high  values 
that  the  system  becomes  unstable. ^ He  further  divided  the  Instability 
mechanism  Into  two  different  physical  forms:  (1)  local  volume  burning 
at  the  end  of  the  cavity,  and  (2)  local  destruction  of  the  propellant. 

In  his  theoretical  analysis,  steady-state  conditions  were  assumed,  but 
this  is  not  practical  in  a real  case.  The  ignition  criterion  was  not 
established  in  the  analysis,  and  the  continuity  equation  did  not  include 
the  mass  addition  due  to  burning.  Friction  loss  in  the  momentum  equa- 
tion was  also  neglected.  Energy  release  due  to  combustion  and  energy 
loss  from  the  hot  gas  to  the  unburned  propellant  was  ignored.  Thus,  the 
analysis  cannot  predict  the  flame  propagation  and  the  actual  combustion 
phenomena  Inside  a propellant  crack. 

Margolin,  et  al. , indicated  that  the  spontaneous  penetration  of 
combustion  gases  into  the  pores  of  a gas-permeable  explosive  charge 
occurs  when  the  Andreev  number  (An)  is  beyond  its  critical  value.®  The 
boundary  conditions  were  found  to  be  important  in  flame  propagation. 
Godal,  in  his  experiment.  Indicated  that  a threshold  diameter  or  a 
critical  width  of  a uniform  cavity  exists,  below  which  flame  will  not 
propagate  into  the  crack.®. 

Krasnov,  et  al. , investigated  the  rate  of  penetration  of  combus- 
tion into  the  pores  of  an  explosive  charge.^®  The  experiment  was  per- 
formed in  a constant-pressure  bomb.  The  experimental  results  showed 
that  the  flame  spreading  rate  is  not  equal  to  the  rate  of  penetration  of 
the  combustion  products  into  the  channel.  Their  theoretical  model  was 
based  on  the  following  crude  assumptions:  (1)  constant  thermal  wave 
penetration  depth  in  the  heated  propellant  along  the  flow  direction. 


^G.  P.  Cherepanov.  "Combustion  in  Narrow  Cavities,"  J APPL  MECH  & 
TECH  PHY,  Vol.  11,  No.  2,  1970,  pp.  276-281. 

®A.  D.  Margolin  and  V.  M.  Margulls.  "Penetration  of  Combustion 
into  an  Isolated  Pore  in  an  Explosive,"  Combustion,  Explosion  & Shook 
Vaves,  Vol.  5,  No.  1,  1969,  pp.  15-16. 

®T.  Godal.  "Flame  Propagation  into  the  Crack  of  a Solid -Propellant 
Grain,"  AMER  INST  AERONAUT  & ASTRONAUT  J,  Vol.  8,  No.  7,  July  1970,  pp. 
1322-1327. 

^®Yu.  K.  Krasnov,  et  al.  "Rate  of  Penetration  of  Combustion  into 
the  Pores  of  an  Explosive  Charge,"  Combustion,  Explosion  & Shock  Waves, 
Vol.  6,  No.  3,  1970,  pp.  262-265. 
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(2)  uniform  temperature  profile  inside  the  heated  propellant  layer,  and 

(3)  steady-state  assumptions  for  the  heat  balance  equation.  Hence, 
their  theoretical  approach  is  Inadequate  for  predicting  flame  spreading. 

Jacobs,  et  al.,  studied  the  pressure  distributions  in  burning 
cracks  which  simulate  the  separation  or  debond  of  solid  propellants  from 
the  motor  casing. Wedge-shaped  channels  with  debond  angles  of  4.75, 
3.17,  and  1.90  degrees  were  tested  in  a combustion  chamber  at  various 
levels  of  pre-pressurization.  Pressure  measurements  were  obtained  at 
three  locations  along  the  debond  and  at  one  location  in  the  pressure 
chamber.  The  propellant  was  doped  with  a highly  flammable  mixture  of 
titanium  and  ammonium  perchlorate  powder  so  that  ignition  under  in- 
creased pressures  could  be  safely  assumed  to  approach  that  of  instan- 
taneous ignition  over  the  entire  surface.  This  ignition  procedure  was 
perhaps  advantageous  for  them  to  compare  the  experimental  results  with  a 
quasi-steady  theoretical  model.  The  actual  ignition  and  flame  spreading 
processes  were  avoided.  Therefore,  their  analytical  results  cannot  be 
used  to  predict  any  transient  behavior  in  crack  combustion. 

Other  studies  of  combustion  inside  propellant  cracks  can  be  found 
in  the  literature. Although  many  interesting  experimental  results 
have  been  obtained,  no  sound  theoretical  model  was  developed  to  describe 
the  transient  processes  governing  crack  combustion. 


^ ^Air  Force  Rocket  Propulsion  Laboratory.  An  Expevimentat  Study  of 
the  Pressure  Distribution  in  Burning  Flaws  in  Solid  Propellant  Grains, 
by  H.  R.  Jacobs,  et  al.,  University  of  Utah.  Edwards,  Calif.,  AFRPL, 
October  1972.  (AFRPL-TR-72-108,  publication  UNCLASSIFIED.) 

^^A.  F.  Belyaev,  et  al.  "Consequences  of  the  Penetration  of  Combus- 
tion into  an  Individual  Pore,"  Combustion,  Explosion  & Shook  Haves,  Vol. 
6,  No.  2,  April-June  1970,  pp.  149-153. 

^^E.  C.  Francis,  G.  H.  Lindsey,  and  R.  R.  Parmerter.  "Pressurized 
Crack  Behavior  in  Two-Dimensional  Rocket  Motor  Geometries,"  AMIP.  INST 
AERONAUT  ASTRONAUT  J,  Vol.  9,  No.  6,  June  1972,  pp.  415-419. 

V.  Kirsanova  and  0.  I.  Leipunskii.  "Investigation  of  the 
Mechanical  Stability  of  Burning  Cracks  in  a Propellant,"  Combustion, 
Explosion  & Shook  Haves,  Vol.  6,  No.  1,  1970,  pp.  68-75. 
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■1 

1 THEORETICAL  MODELING 


The  physical  model  chosen  to  simulate  a crack  in  a propellant  grain 
is  depicted  schematically  in  Figure  1.  It  represents  a pore  located  in 
the  transverse  direction  to  the  main  flow  of  the  rocket  chamber.  The 
crack  can  be  of  variable  geometry  with  the  opening  initially  starting  at 
or  near  the  propellant  surface  which  is  exposed  to  the  main  chamber 
flow.  This  physical  model  can  also  be  easily  extended  to  simulate 
conditions  of  debond  or  separation  of  the  propellant  from  the  motor 
casing. 


FIGURE  1.  The  Physical  Model  of  a Transverse 
Crack  in  Solid  Propellant. 

The  following  basic  assumptions  are  used  to  facilitate  a theoreti- 
cal description  of  the  transient  combustion  processes  inside  a crack: 

1.  All  chemical  reactions  occur  near  the  propellant  crack  surface, 
in  a combustion  zone  which  is  so  thin  that  it  can  be  considered  as  a 
plane.  That  is,  the  thickness  of  the  combustion  zone  above  the  propel- 
lant surface  is  much  smaller  than  the  gap  width  in  the  propellant  crack. 
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2.  Rate  processes  at  the  propellant  surface  are  quasi-steady  in 
the  sense  that  the  characteristic  times  associated  with  the  gaseous 
flame  and  preheated  propellant  are  short  compared  to  that  of  the  pres- 
sure transient  variation. 

3.  The  gases  flowing  in  the  propellant  crack  obey  the  Clausius  or 
Noble-Abel  equation  of  state.  This  dense  gas  relation  can  adequately 
describe  the  departure  from  the  ideal  gas  law  at  high  pressures. 

4.  The  bulk  flow  in  the  pore  is  one-dimensional.  Recent  experi- 
ments^^ conducted  by  Kuo,  et  al. , have  revealed  that  the  one-dimensional 
assumption  is  appropriate  for  flow  in  cracks  when  the  pressure  gradient 
along  the  crack  is  greater  than  a critical  value  which  normally  is 
extremely  small. 

5.  The  compressibility  of  the  solid  propellant  is  neglected.  The 
change  of  port  area  in  propellant  cracks  is  solely  due  to  the  propellant 
surface  regression.  The  authors  recognize  that  the  assumption  of  negli- 
gible propellant  deformation  may  be  unrealistic.  However,  before  this 
effect  can  be  adequately  treated,  a reliable  transient  combustion  model 
has  to  be  developed  to  supply  the  necessary  temporal  and  spatial  pres- 
sure distributions  required  by  an  appropriate  stress-analysis  code. 


GOVERNING  EQUATIONS  FOR  THE  GAS  PHASE 

To  describe  the  transient  gas-phase  behavior  inside  a propellant 
crack,  we  derive  the  unsteady,  quasi-one-dlmenslonal  mass,  momentum  and 
energy  equations  from  a balance  of  fluxes  which  act  on  an  elemental 
control  volume  as  shown  in  Figure  1.  The  resultant  conservation  equa- 
tions expressed  in  divergence  form  are: 

Mass  Conservation 


9(pA  ) 


3(p  u A ) 

EL  _ 


'b  '’pt  ^b 


(1) 


K.  Kuo,  A.  T.  Chen,  and  T.  R.  Davis.  "Transient  Flame  Spread- 
ing and  Combustion  Processes  inside  a Solid  Propellant  Crack,"  presented 
at  15th  AIAA  Aerospace  Soienaes  Meeting,  January  1977.  (AIAA  Paper  No. 
77-14) 
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Momentum  Conservation 
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Energy  Conservation 
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(3) 


The  conservation  equations  are  simplified  further  by  an  order-of- 
magnltude  analysis  which  shows  that  the  following  terms  are  negligible: 

1.  Forces  between  gas  molecules  due  to  the  viscous  normal  stress 
In  the  axial  direction,  created  by  the  axial  velocity  gradient 
In  the  momentum  equation 

2.  Viscous  dissipation  and  rate  of  work  done  by  the  force  due  to 
viscous  normal  stresses  In  the  energy  equation 

3.  Axial  heat  conduction  between  gas  molecules  In  the  energy 
equation. 

After  rearranging  the  conservation  equations,  a set  of  three  govern- 
ing equations  Is  obtained,  l.e.,  velocity-variation,  temperature- 
variation,  and  pressure-variation  equations  which  are  shown  below. 
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Velocity-Variation  Equation 
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Temperature-Variation  Equation 
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Pressure-Variation  Equation 
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The  three  governing  equations  are  first  order,  coupled.  Inhomoge- 
neous, and  non-llnear  partial  differential  equations.  The  unknowns 
which  vary  strongly  with  respect  to  the  spatial  variables  are:  gas 
temperature,  pressure,  velocity  and  the  temperatures  at  the  propellant 
surface  and  Inert  wall  If  present.  There  are  other  unknown  parameters 
which  vary  weakly  with  respect  to  spatial  variables.  These  parameters 
Include  the  port  area  of  the  crack,  the  burning  rate  of  the  propellant, 
the  convective  heat  transfer  and  friction  coefficients,  the  gas  flow 
angle,  the  local  blowing  velocity  normal  to  the  propellant  surface,  the 
wetted  perimeter  and  the  burning  perimeter  at  various  locations  along  the 
crack. 


In  order  to  solve  all  the  unknowns  mentioned  above,  additional 
equations  and  boundary  conditions  together  with  some  physical  Input 
functions  must  be  specified.  The  equation  which  governs  the  surface 
temperature  variation  of  the  solid  phase  Is  described  below. 


EQUATION  FOR  DETERMINATION  OF  PROPELLANT  SURFACE  TEMPERATURE 

The  propellant  surface  temperature  at  a fixed  location  along  the 
crack  before  the  attainment  of  Ignition  Is  calculated  from  the  solid- 
phase  heat  conduction  equation  which  Is  written  In  the  unsteady,  one- 
dlmenslonal  form 


3 T (t,  y) 


3 T (t,  y) 


where  the  length  variable,  y.  Is  measured  perpendicular  to  the  local 
propellant  crack  surface. 

Since  this  equation  Is  first  order  In  time  derivative  and  second 
order  In  space  derivative,  the  following  Initial  and  boundary  conditions 
are  required 

V ' %i 
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An  approximate  solution  of  the  heat  equation,  which  Is  quite  accurate 
and  saves  much  computation  time  may  be  obtained  by  using  the  Integral 
method. Applying  this  method,  a quantity  (called  the  thermal 

penetration  distance)  is  defined,  so  that  the  y > 6^,  the  solid  is  at 
equilibrium  temperature,  and  there  Is  no  heat  transfer.  It  Is  illus- 
trated In  Figure  2. 


FIGURE  2.  Temperature  Profile  In  the  Solid  Propellant. 


A third  degree  polynomial  is  used  with  time  varying  coefficients  to 
approximate  the  transient  temperature  profile  in  the  solid  phase 

T = k^(t)  + k (t)  y + k (t)  y^  + k (t)  y^  (11) 

pr  u 1 l j 


^^T.  R.  Goodman.  "The  Heating  of  Slabs  with  Arbitrary  Heat  Inputs," 
J AEROSPACE  SCI,  Vol.  26,  March  1959,  pp.  183-188. 

^^T.  R.  Goodman.  "Application  of  Integral  Methods  to  Transient 
Non-linear  Heat  Transfer,"  Advanoee  in  Heat  Transfer,  Vol.  1,  1964, 
Academic  Press,  New  York,  pp.  51-122. 
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where  the  coefficients  Icq,  k2,  and  k3  may  be  time^dependent.  To 
determine  them,  four  boundary  conditions  are  necessary.  These  are 


T 

pr 


(t,  6^)  - T 


pi 


(12) 
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(t,  6^) 
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(13) 
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(t,  0)  = - ^ [T(t)  - 
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ps 


(14) 


(15) 


After  applying  the  integral  method,  the  resulting  ordinary  differ- 
ential equation  describing  the  rate  of  increase  of  propellant  surface 
temperature  is 


d T 

E£  = 

dt 


4 o h ^(T  - T 
pr  c gs 

(T  - T J(2T  - T 
ps  pi 


ps 


- V 


(16) 


The  initial  condition  used  for  this  equation  is:  T 
is  a small  value  used  to  avoid  singularity.  ^ 


T + e,  where  e 
pi 


EQUATION  OF  STATE  FOR  GAS  AND  SOLID  PHASES 

In  addition  to  the  governing  equations  in  the  gas  phase  and  the 
equation  describing  the  rate  of  Increase  in  Tpg,  the  equation  of  state 
for  the  gas  and  solid  phases  must  be  specified.  As  discussed  earlier, 
the  Clausius  or  Noble-Abel  equation  of  state 

P(l/p  - b)  - RT  (17) 

was  chosen  to  account  for  any  non-ideal  gas  behavior.  The  term,  b, 
represents  the  co-volume  or  volume  occupied  by  the  individual  gas  mole- 
cules. One  obvious  advantage  of  using  the  Noble-Abel  equation  of  state 
is  its  flexibility;  it  can  be  readily  reduced  to  the  ideal  gas  equation 
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by  setting  the  value  of  the  co-volume  to  zero.  The  statement  of  con- 
stant density  for  the  solid  propellant  serves  as  the  equation  of  state 
for  the  solid  phase. 


EMPIRICAL  HEAT  TRANSFER  AND  FRICTION  COEFFICIENT  CORRELATION 


_ The  expression  for  the  local  convective  heat-transfer  coefficient, 
h^,  is  deduced  from  the  conventional  Dittus-Boelter  correlation  for 
turbulent  flow  in  pipes. Variation-of  the  physical  properties  of  the 
gas  across  the  boundary  layer  is  considered  by  evaluatlng_the  properties 
at  an  average  film  temperature,  T^f.  The  expression  for  h^,,  used  in  the 
model,  is 


h = 0.0346  Pr”®‘^  c (P  u/R)°'®  V T (x  d (18) 

c p k af  h 


Before  local  ignition  h = h = h , whereas  after  local  ignition 
— — — ^ CD  cw  r 

h =0  and  h = h . 
cp  cw  c 


The  above  equation  has  incorporated  the  temperature  dependence  of 
the  viscosity  which  was  determined  from  the  calculations  of  Svehla  for 
the  viscosity  of  air  at  high  temperature. ^ ^ It  was  shown  by  Bartz  that 
Svehla 's  values  are  very  close  to  those  in  the  previously  published  NBS 
data.^*^’^^  The  coefficient  of  viscosity-temperature  relation,  is 
defined  as 


V,  = 8.699  X 10“^  M 
k w 


(19) 


^®A.  Peretz,  et  al.  "Starting  Transient  of  Solid  Propellant 
Rocket  Motors  with  High  Internal  Gas  Velocities,"  AMER  INST  AERONAUT 
ASTRONAUT  J,  Vol.  11,  No.  12,  December  1973,  pp.  1719-1727. 

'^National  Aeronautics  and  Space  Administration.  Estimated  Vis- 
aosities  and  Thermal  Conductivities  of  Gases  at  High  Temperatures,  by  R. 
A.  Svehla.  Cleveland,  Ohio,  NASA,  1962.  (Publication  UNCLASSIFIED.) 

^®D.  R.  Bartz.  "Survey  of  the  Relationship  Between  Theory  and 
Experiment  for  Convective  Heat  Transfer  from  Rocket  Combustion  Gases," 
Advances  in  Tactical  Rocket  Propulsion,  AGARD  Conference  Proceedings  No. 
1.  Technivision  Services,  Maidenhead,  England,  August  1968,  pp.  291- 
381. 

Hilsenrath,  et  al.  Tables  of  Thermal  Properties  of  Gases, 

U.S.  Department  of  Commerce,  NBS  Circular  564,  November  1955. 
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The  Prandtl  number  is  calculated  from  Svehla's  equation^ 


1,77  Y - 0.45 


The  correlation  for  the  friction  coefficient  used  in  this  study  is 


0.4491  (d  /x) 
n 


E /d.  1. 

3.7 


nTT  a.  1*65  , .0.05' 

46  RV  T (d. /x) 

K n 

.1  0-5 

p u dh  Cf 


where  eg/d},  is  the  relative  equivalent  sand  roughness.  This  expression 
is  actually  a modified  form  of  the  well-known  Colebrook  formula^^*^^ 
for  turbulent  flow  in  pipes  with  roughness.  Entrance  effects  are  taken 
into  account  by  modifying  the  friction  coefficient  by  a power  function 
of  the  distance-to-dlameter  ratio;  the  original  Colebrook  expression  is 
obtained  for  s/d^  > 20.  Additionally,  Eq.  (21)  implicitly  accounts  for 
the  temperature  dependence  of  the  gas  density  and  viscosity,  each  evalu- 
ated at  the  average  film  temperature,  T . 

ai 

After  the  propellant  surface  starts  burning  locally,  the  value  of 
the  friction  coefficient  is  set  to  zero  due  to  the  attenuation  of  wall 
shear  stress  caused  by  surface  blowing. 


PROPELLANT  BURNING  RATE  EXPRESSION 

In  the  absence  of  a more  adequate  erosive  burning  rate  expression, 
the  Lenoir-Robillard  burning  rate  formula^**  is  used  temporarily  to 
account  for  the  erosive  burning  effect 

r,  = a p”  + K h exp(-6  r.  p /u  p)  (22) 

b e c b pr 

Schlicting.  "Turbulent  Flow  Through  Pipes,"  Boundary  Layer 
Theory,  6th  Ed.,  McGraw-Hill,  New  York,  1968. 

F.  Colebrook.  "Turbulent  Flow  in  Pipes  with  Particular  Ref- 
erence to  the  Transition  Region  Between  the  Smooth  and  Rough  Pipe  Laws," 
J Inst  Civil  Engineers,  Vol.  11,  pp.  133-156,  1938-39. 

^“*3.  M.  Lenoir  and  G.  Roblllard.  "A  Mathematical  Method  to  Pre- 
dict the  Effects  of  Erosive  Burning  in  Solid  Propellant  Rockets,"  6th 
Symp.  (Int.)  on  Combustion,  The  Combustion  Institute  (1956),  pp.  663-7. 
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This  relation  Is  actually  a modified  form  o^  the  Lenolr-Robillard  formula, 
since  the  original  expression  proposed  for  h(.  Is  based  on  the  Chllton- 
Colburn  correlation  for  turbulent  flow  over  a flat  plate,  whereas  the 
expression  for  h^.  In  Eq.  (22)  Is  the  modified  Dlttus-Boelter  correlation 
(Eq.  (18)).  The  local  burning  rate  as  given  by  the  above  relation  Is 
strongly  coupled  to  the  gas  dynamics  inside  the  crack. 

PROPELLANT  IGNITION  CRITERION 

After  Tps  Is  solved,  a simplified  Ignition  temperature  criterion  Is 
used  to  determine  the  burning  condition  of  the  solid  propellant  along 
the  crack.  In  order  to  avoid  the  step  function  change  In  burning  rate, 
a two-temperature  criterion  Is  used  to  achieve  full  Ignition  within  a 
finite  time  Interval;  when  Tpg  Is  equal  to  the  propellant  starts 

to  gasify.  As  soon  as  Tpg  reaches  l^gn*  full  Ignition  condition  is 
reached. 


INITIAL  AND  BOUNDARY  CONDITIONS 

The  initial  conditions  required  for  the  solution  of  the  system  of 
governing  equations  are; 

u(0,  x)  = u^ 

T(0,  x)  = T. 

P(0.  x)  = P^  (23) 

where  Uj^,  Tj^,  and  Pj^  are  the  initial  velocity,  temperature,  and  pressure 
In  the  propellant  crack. 

The  number  of  boundary  conditions  required  depends  upon  the  flow 
direction  and  Mach  number  at  the  opening  of  the  crack.  When  the  gas  In 
the  rocket  chamber  flows  subsonlcally  Into  the  propellant  crack,  two 
boundary  conditions  at  the  opening  of  the  propellant  crack  are  specified 
to  represent  the^  Interface  conditions  at  the  crack  entrance;  and  one 
boundary  condition  at  the  closed  end  of  the  propellant  crack  Is  speci- 
fied to  Indicate  that  the  gas  velocity  Is  zero.  In  the  mathematical 
form  they  are  represented  as  follows: 


P(t,  0)  - P (t) 
c 

T(t,  0)  - T (t) 
c 

U(t,  Xj^)  * 0 (24) 
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When  the  flow  of  Che  gas  out  of  Che  crack  Is  subsonic,  two  boundary 
conditions  are  specified.  They  are  written  as 

P(t,  0)  - P (t) 
c 

u(t,  - 0 (25) 

When  the  outflowing  gas  Is  supersonic,  only  one  boundary  condition  can 
be  specified,  that  Is 

u(t,  - 0 (26) 

In  order  Co  solve  the  system  of  governing  equations  with  a second 
order  numerical  scheme,  extraneous  boundary  conditions  are  required  In 
addition  to  the  physical  ones.  The  procedure  to  attain  these  conditions 
Is  discussed  In  the  following  section. 

NUMERICAL  SOLUTION  TECHNIQUE 

The  governing  equations  (Eq.  (4)  through  (6)),  describing  the 
dynamic  behavior  of  a Noble-Abel  dense  gas  In  a propellant  crack  were 
found  to  be  totally  hyperbolic. The  numerical  techniques  developed 
were  chosen  on  the  basis  of  stability,  maximum  accuracy,  and  computa- 
tional efficiency.  Recently  obtained  experience  in  solving  hyperbolic 
partial  differential  equatlons^^  was  utilized,  and  a generalized  Impli- 
cit scheme^^*^®  based  on  central  differences  In  space-wise  derivatives 
was  chosen  to  solve  the  governing  equations.  A quasl-llnearlzatlon 
method,  used  to  linearize  the  Inhomogeneous  terms  in  the  governing 
equation,  was  combined  with  a stable  predictor-corrector  procedure.^® 
Finally,  the  governing  equations,  in  their  finite-difference  form,  were 
arranged  in  a block-tridiagonal  matrix  form,^®  which  allowed  an  efficient 
computation. 

^®R.  Courant  and  D.  Hilbert.,  Methods  of  Mathematical  Physics ^ Vol. 
2,  Interscience  Publishers,  Inc.,  New  York,  July  1966,  pp.  407-550. 

K.  Kuo.  "Theory  of  Flame  Front  Propagation  In  Porous  Propel- 
lant Charges  under  Confinement,"  Ph.D.  Thesis.  Princeton  University, 
Princeton,  N.  J. , August  1971.  (AMS  Report  T-1000. ) 

^^R.  D.  Rlchtmyer  and  D.  w.  Morton.  Difference  Methods  for  Initial- 
Value  Problems.  Interscience  Publishers,  New  York,  1967. 

^®M.  G.  Salvadorl  and  M.  L.  Baron.  Numerical  Methods  in  Engineer- 
ing. Englewood  Cliffs,  N.  J.,  Prentice  Hall,  1961. 

^®A.  T.  Chen.  "Theoretical  Modeling  and  Numerical  Solution  of 
Transient  Flame-Spreading  in  Solid-Propellant  Cracks,"  M.S.  Thesis, 
Mechanical  Engineering  Department,  Pennsylvania  State  University,  May 
1976. 
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The  use  of  central  difference  formulation  for  all  spacewlse  deriva- 
tives in  the  governing  equations  requires  six  boundary  conditions  for 
the  solution.  This  implies  that  when  the  gas  flows  into  the  cracic, 
three  extraneous  boundary  conditions  are  required  in  addition  to  the 
three  physical  ones  described  by  Eq.  (24).  The  three  extraneous  boundary 
conditions  are  used  to  determine  the  flow  velocity  at  the  opening  of  the 
crack,  and  the  pressure  and  temperature  at  the  closed  end  of  the  crack. 

When  the  gas  flows  out  of  the  crack  subsonlcally,  in  addition  to 
the  two  physical  boundary  conditions  described  by  Eq.  (25),  two  extran- 
eous boundary  conditions  are  required  at  the  opening  of  the  crack  for 
the  determination  of  the  pressure  and  temperature.  Also,  two  extraneous 
boundary  conditions  are  required  at  the  closed  end  of  the  crack. 

When  the  outflowing  gas  is  supersonic,  in  addition  to  the  physical 
boundary  condition  given  by  Eq.  (26),  three  extraneous  boundary  condi- 
tions are  needed  at  the  opening  of  the  crack  and  two  other  extraneous 
boundary  conditions  are  needed  at  the  tip  of  the  crack. 

These  extraneous  boundary  conditions,  which  are  also  frequently 
called  compatibility  relations,  are  given  below: 


(27) 


(28) 


s 


where  F^,  Fj,  and  Fp  represent  th6  inhomogeneous  terms  on  the  right  hand  j 

side  of  Eq.  (4),  (5),  and  (6).  The  subscripts  I,  II,  and  III  designate  j 

the  right-running,  left-running  characteristic  lines  and  the  particle  ! 

path  line,  respectively.  These  relations  were  obtained  by  transforming  j 

the  governing  hyperbolic  equations  into  their  characteristic  form.  j 

Together  with  the  physical  boundary  conditions,  the  compatibility  rela-  | 

tions  form  a closed  system  for  the  determination  of  velocity,  tempera- 
ture, and  pressure  at  both  ends  of  the  crack.  ’ 

The  extraneous  boundary  conditions  are  numerically  evaluated  by  a 
fourth  order  Runge-Kutta  integration  technique.  In  previous  calcula- 
tions, the  integration  of  each  compatibility  relation  was  treated  con- 
secutively. Currently,  these  equations  are  simultaneously  Integrated  to 
Improve  the  numerical  stability  of  the  solution. 
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EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


TEST  RIG  DESIGN 

An  overall  view  of  the  experimental  combustion  chamber  used  to 
study  the  flame  spreading  and  pressurization  phenomena  in  propellant 
cracks  is  shown  in  Figure  3.  The  design  of  this  test  rig  was  guided  by 
the  following  primary  considerations: 


1 ■ 

[ FIGURE  3.  Overall  View  of  the  Assembled  Gaseous  Ignition 

System  Mounted  With  the  Combustion  Chamber  on  the  Test 
t , Stand. 


21 


T 


- 


11.^11 1 


NWC  TP  5943 

1.  The  experimental  set-up  should  be  as  close  as  possible  to  the 
physical  conditions  considered  in  the  analytical  model  so  that  the 
experimental  results  can  be  used  directly  for  the  verification  of  the 
theoretical  solutions. 

2.  The  test  rig  should  have  sufficient  flexibility  to  permit  a 
vide  range  of  variation  in  propellant  crack  geometry,  properties,  chamber 
pressure,  pressurization  rate,  igniter  gas  temperature  and  the  main 
chamber  mass  flow  rate  in  the  transverse  direction  to  the  crack.  This 
flexibility  in  testing  is  essential  for  the  establishment  of  a funda- 
mental data  base  in  propellant  crack  combustion. 

A schematic  diagram  of  the  test  rig  is  shown  in  Figure  4,  and  an 
assembly  view  of  the  propellant  crack  combustion  chamber  is  shown  in 
Figure  5.  The  propellant  containing  a simulated  crack  is  cast  into  a 
brass  mold.  A special  mandrel  is  used  during  the  propellant  molding 
process  to  produce  a 2.54  cm  deep,  channel-shaped  crack  of  length,  L, 
and  gap  width,  6.  Thus  far,  two  different  crack  geometries  have  been 
used:  for  the  large  crack,  L = 20  cm  and  the  average,  6 = 0.52  cm;  for 
the  small  crack,  L = 15  cm  and  the  average,  6 = 0.08  cm.  The  main  olock 
of  the  combustion  chamber  has  five  parts  machined  into  its  backside 
along  the  longitudinal  axis  of  the  crack  to  accept  four  piezoelectric 
pressure  transducers  and  a safety  burst  diaphragm.  The  pressure  trans- 
ducers are  located  at  a distance  of  0.0,  4.8,  13.8,  and  18.8  cm  from  the 
initial  crack  opening.  The  window  assembly  consists  of  a 1/4-inch 
sacrificial  Plexiglas  plate  and  a main  Plexiglas  window.  This  window 
assembly  along  with  a stainless  steel  retainer  is  used  to  close  the 
combustion  chamber  and  permit  visual  recording  of  the  luminous  flame 
front  inside  the  propellant  crack.  An  interchangeable  exit  nozzle  is 
connected  to  the  combustion  chamber  for  control  of  the  rate  of  chamber 
pressurization.  The  combustion  chamber  was  designed  for  a maximum 
static  pressure  of  800  atmospheres. 

An  ignition  system  has  been  developed  using  spark  plugs  to  ignite  a 
gaseous  mixture  of  oxygen  and  methane.  A schematic  drawing  of  the 
igniter  system  is  shown  in  Figure  6.  Both  fuel  and  oxidizer  feedlines 
contain  a pressure  regulator  and  critical  flow  orifice  which  permit 
exact  control  of  the  fuel-oxidizer  flow  rates.  The  gaseous  reactants 
are  first  mixed  in  the  igniter  chamber  by  tangentially  injecting  methane 
into  the  oxygen  flow  stream.  The  mixed  reactants  then  enter  the  test 
section  through  a rectangularly  oriented,  multi-perforated,  convergent 
nozzle.  This  versatile  igniter  design  provides  the  necessary  varia- 
bility to  produce  a wide  range  of  igniter  strength. 
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(A)  Main  Chamber  Block  (D)  Main  Window 

(B)  Brass  Mold  Containing  (E)  Window  Retainer 

Propellant-Crack  Specimen 

(C)  Sacrificial  Window 


FIGURE  5.  Assembly  View  of  the  Combustion  Chamber. 


IGNITION  CHAMBER 


SPARK  PLUGS 


FIGURE  6.  Schematic  Diagram  of  the  Gaseous  Ignition  System. 
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INSTRUMENTATION 


A block  diagram  of  the  data  acquisition  system  is  shown  in  Figure  7. 
The  components  of  this  instrumentation  system  are: 

1.  A high-speed  magnetic  tape  recorder  (Hewlett  Packard  model 
3924B)  with  a simultaneous  14  channel  record/reproduce  capability  at  a 
rate  up  to  60  in/s.  Seven  channels  record  in  FM  mode,  and  the  other 
seven  channels  record  in  the  direct  or  AM  mode. 

2.  An  analog-to-dlgital  converter  (Blomation  1015  Transient  Wave- 
form Digitizer)  provides  4-channel  data  storage  with  1024  word  memory 
per  channel,  and  is  capable  of  a real-time  resolution  down  to  10  ys. 

The  output  of  the  analog-to-digital  converter  is  used  to  control  the 
horizontal  and  vertical  deflection  of  an  X-Y  plotter  to  obtain  a hard 
copy  of  the  transient  data. 

3.  A 16  mm  high-speed  camera  (FASTAX  Model  W-163269)  is  used  with 
a 100-foot  film  roll  to  record  the  propagation  of  the  luminous  flame 
front  and  regression  of  the  propellant  surface. 


FIGURE  7.  Block  Diagram  of  the  Instrumentation  System. 
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4.  Piezoelectric  pressure  transducers  (Kistler  Model  601B) 
equipped  with  water-cooled  adapters  (Kistler  Model  628C)  are  used  to 
measure  the  pressure  at  four  locations  along  the  propellant  crack. 

Charge  amplifiers  (Kistler  Model  504  E Dual  Mode  Amplifiers)  are  used  to 
condition  the  pressure  signal  for  tape  recording.  These  transducers  are 
recess-mounted  in  the  combustion  chamber  and  can  be  used  to  measure 
pressures  up  to  1000  atmospheres. 


TEST  PROCEDURE 

A typical  test  firing  is  conducted  in  the  following  manner.  In 
order  to  prevent  unburned  Igniter  gases  from  entering  the  crack,  a strip 
of  tape  is  used  to  cover  the  entrance  of  the  crack.  The  fuel  and  oxi- 
dizer flow  rates  are  set  to  predetermined  values.  As  the  cold  igniter 
gases  pass  over  the  sealed  crack  entrance,  a high-speed  camera  is  turned 
on.  A timing  signal  is  then  activated  which  places  light  pulses  on  the 
film  and  simultaneously  sends  a pulse  to  the  tape  recorder  to  synchro- 
nize the  film  with  the  pressure  data  on  the  tape.  Ignition  is  then 
achieved  with  two  conventional  automobile  spark  plugs.  Hot  igniter- 
gases  pass  from  the  igniter  chamber  through  a multi-perforated  nozzle  to 
the  combustion  chamber;  the  tape  covering  the  crack  entrance  burns  away 
immediately  and  the  hot  gases  enter  the  crack. 

After  an  induction  interval,  the  propellant  at  the  crack  entrance 
reaches  ignition  condition.  The  flame  spreading  process  is  then  observed 
by  the  16  mm  high-speed  camera.  The  pressurization  process  in  the 
propellant  crack  is  simultaneously  recorded  on  the  magnetic  tape. 

DISCUSSION  OF  RESULTS 


Calculations  have  been  made  with  the  present  theoretical  model  to 
compare  with  initial  experimental  data  obtained  from  the  crack  combus- 
tion facility.  Additionally,  a parametric  study  was  conducted  with  the 
theoretical  model  to  ascertain  the  Influence  of  variations  in  chamber 
pressurlzaton  rate,  crack  gap  width,  and  ignition  parameters  on  the 
flame  spreading  and  resultant  pressurization  Inside  propellant  cracks. 


EXPERIMENTAL  DATA  AND  THEORETICAL  COMPARISON 

A limited  number  of  firings  with  a smokeless  propellant  have  been 
conducted  in  the  newly  developed  crack  combustion  facility.  The  propel- 
lant used  in  these  firings  was  furnished  by  NWC.  As  described  earlier 
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in  this  report,  two  different  crack  geometries  have  been  tested.  The 
larger  crack  (6  = 0.5  cm)  was  used  during  the  phase  of  initial  system 
checkout  and  will  be  retested  for  documentation  in  a later  report.  The 
experimental  results  of  a smaller  crack,  having  a length  of  15  cm  and  a 
gap  width  of  0.051  cm  at  the  entrance  with  a gradual  increase  to  0.107 
cm  at  the  tip,  are  discussed  herein.  (Although  we  call  the  above  con- 
figuration a "small"  crack,  the  gap  width  may  be  an  order  of  magnitude 
larger  than  that  of  the  actual  cracks  found  in  rocket  grains.) 

Typical  pressure-time  traces  recorded  from  a test  firing  of  a 
"small"  crack  configuration  are  shown  in  Figure  8.  G2  represents  the 
pressure  history  at  the  entrance  of  the  crack.  G5  was  initially  covered 
by  the  propellant  specimen  and  therefore  senses  pressure  at  a later 
time.  In  this  test  a large  exhaust  nozzle  was  used  to  produce  a rela- 
tively low  pressurization  rate  at  the  crack  entrance.  This  resultant 
pressurization  rate  is  close  to  the  typical  value  obtained  during  the 
starting  transient  of  many  rocket  motors.  An  examination  of  the  G2 
through  G4  pressure-time  traces  reveal  that  the  spatial  pressure  gradi- 
ent inside  the  crack  is  very  low.  This  nearly  uniform  spatial  pressure 
history  is  due  to  the  relatively  low  value  of  dPc/dt  and  to  a signifi- 
cantly large  gap  width;  these  conditions  are  favorable  for  the  reduction 
of  pressure  wave  phenomena.  As  shown  in  Figures  9 and  10,  the  predicted 
P-t  profiles  at  two  stations  inside  the  crack  are  in  close  agreement 
with  the  experimental  data. 


0 0.8  1.0  1.8  2.0 

TIME.  ■ 

FIGURE  8.  Experimental  Pressure-Time  Traces  from  DNC  Test  No.  11. 
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FIGURE  9.  Comparison  of  Theoretical  and  Experimental  P-t 
Profiles  at  x = 4.8  cm  (Simulation  of  DNC  Test  No.  11). 


TIME,  t 

FIGURE  10.  Comparison  of  Theoretical  and  Experimental  P-t 
Profiles  at  X ■ 13.8  cm  (Simulation  of  DNC  Test  No.  11). 
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In  this  test  an  exceedingly  long  Induction  period  was  required  for 
the  Igniter  gases  to  burn  away  the  protective  tape  and  establish  Igni- 
tion at  the  crack  entrance.  Consequently,  no  useful  flame  spreading 
data  was  obtained.  However,  In  another  firing  with  nearly  Identical 
test  conditions,  the  luminous  flame  front  was  found  to  propagate  slowly 
near  the  crack  entrance  and  then  suddenly  accelerate  to  a nearly  uniform 
velocity.  The  film  showing  the  characteristic  features  of  the  flame 
spreading  process  was  presented  at  the  last  High  Energy  Propellant 
Safety  (HEPS)  Meeting  held  at  NSWC/WOL,  December  1976.  A similar  be- 
havior In  the  flame  spreading  rate  was. also  observed  experimentally  by 
Balyaev. ^ 

Figure  11  describes  the  calculated  gas  temperature  distribution 
along  the  crack  at  various  times.  Because  of  the  high  flame  temperature 
and  short  transient  times,  no  attempt  was  made  to  measure  the  gas  tem- 
perature Inside  the  crack.  The  Input  data  used  to  make  these  predic- 
tions are  listed  In  Table  1.  The  locus  of  the  calculated  Ignition  front 
Is  superimposed  on  this  figure.  The  gradient  In  the  temperature  distri- 
bution Is  very  pronounced  Indicating  the  penetration  of  hot  gases  Into 
the  propellant  crack.  For  times  greater  than  50  ms,  the  gas  temperature 
at  the  opening  portion  of  the  crack  decreases  slightly;  this  Is  mainly 
due  to  the  outflow  of  product  gases. 


FIGURE  II.  Calculated  Gas  Temperature  Distribution  for 
Various  Times  (Simulation  of  DNC  Test  No.  11) . 
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TABLE  1.  Computer  Program  Input  Variable. 


X 

s 

0.0  cm 

p 

\ 

= 

15.05  cm 

B 

X 

= 

0.0  g /g 

b 

= 

1.0  cm  /g 

M 

V 

= 

23.97  g /g-mole 
r • 

Y 

= 

1.233 

P 
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1.741  g/cm^ 

pr 

-3 

Y 
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0.44  X 10  cal/s-cm-K 

pr 

2 
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0.001  cm  /s 

pr 

e /d 
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0.001 

s h 

^f 
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3000  K 

T . 

780  K 

crl 

T, 

s 

1000  K 

Ign 

'^1 

= 

298  K . 

a 

= 

2.752  X 10  ^ (cm/s)/(g^/cm^)° 

n 

= 

0.683 

K 

— 

15.9  cm^-K/cal 

e 

B 

= 

125 

T 

as 

2800  K 

c 

p 

s 

Experimental  P-t  data  at  G2  (See  Figure  8) 

c 

Figure  12  shows  the  calculated  gas  velocity  profiles  that  exist  In 
the  small  crack  at  various  times.  At  a time  of  1 ms  the  hot  gases  are 
flowing  Into  the  crack  at  a high  velocity  because  of  a large  Initial 
pressure  gradient.  At  25  ms  the  velocity  of  the  Inflowing  gases  de- 
creases as  the  pressure  distribution  along  the  crack  becomes  more  uni- 
form. A negative  velocity  Is  obtained  at  50  ms  In  a region  near  the 
crack  opening  which  Indicates  an  outflow  of  product  gases.  As  the  flame 
front  progresses  along  the  crack,  the  resultant  gasification  processes 


8 0 


DISTANCE  FROM  CRACK  OFENING,  cm 

FIGURE  12.  Calculated  Velocity  Distribution  at 
Various  Times  (Simulation  of  DNC  Test  No.  11). 


cause  an  Increase. in  velocity  of  the  gases  flowing  out  of  the  crack. 

Due  to  the  large  increase  in  port  area  at  opening  of  the  crack  (see 
Figure  13) , the  velocity  of  the  out-flowing  gases  decreases  between  100 
and  200  ms.  During  the  entire  pressurization  process,  the  gas  velocity 
near  the  closed  end  of  the  crack  is  very  low;  this  is  due  to  the  lack  of 
a favorable  pressure  gradient.  This  very  low  flow  velocity  causes  a 
significant  reduction  in  the  convective  heating  of  the  propellant  sur- 
face. Hence,  a much  longer  time  is  required  for  the  propellant  surface 
near  the  crack  tip  to  achieve  a state  of  ignition.  This  delay  is  readily 
evident  in  the  calculated  flame  front  propagation  shown  in  Figure  14. 


PARAMETRIC  STUDY 


A parametric  study  was  conducted  with  the  theoretical  model  to 
determine  the  Influence  of  variations  in  chamber  pressurization  rate, 
crack  gap  width,  and  Ignition  parameters  on  the  flame  spreading  and 
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resultant  pressurization  inside  propellant  cracks.  Unless  otherwise 
described,  the  input  constants  used  in  this  theoretical  study  are  given 
in  Table  1. 

Figure  15  describes  the  pressure  history  at  the  tip  of  a small, 
constant  area  crack  (6  > 0.02  cm)  for  various  chamber  pressurization 
rates.  Three  sets  of  curves  are  shown  in  this  plot  for  dPc/dt  varying 
from  0.5  X 10^  to  1.5  x 10^  atm/s.  Each  set  of  curves  compares  the 
calculated  pressurization  history  at  the  tip  of  the  crack  with  the 
assumed  chamber  pressure  variation.  These  assumed  chamber  pressuriza- 
tion rates  are  of  the  same  order  of  magnitude  that  a crack  would  experi- 
ence during  the  starting  transient  of  a rocket  motor.  As  shown  in  this 
plot,  the  pressurization  rate  at  the  tip  of  the  crack  changes  abruptly 
due  to  the  gasification  of  propellant  surface  inside  the  crack.  The 
time  required  to  attain  this  "rapid  pressurization"  is  aptly  designated 


FIGURE  15.  Influence  of  Rocket  Chamber  Pressurization 
Rate  on  the  Internal  Pressurization  Process  Inside  a 
Narrow  Crack  (6  ■ 0.02  cm). 
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This  time  can  either  be  defined  as:  the  time  when  the  pressure  at 
the  crack  tip  exceeds  the  chamber  pressure,  or;  when  the  two  tangent 
lines,  drawn  on  each  side  of  the  curve  describing  the  rapid  pressure 
transition.  Intersect.  For  small  cracks  and  large  dPc/dt,  both  methods 
for  obtaining  t^^p  are  nearly  equivalent.  As  shown  In  Figure  15,  tgp 
decreases  for  Increasing  chamber  pressurization  rates.  This  trend  Is 
mainly  due  to  the  Increase  of  forced  convection  driven  by  a large  stream- 
wise  pressure  gradient.  This  augmentation  In  the  rate  of  convective 
heat  transfer  significantly  enhances  the  rate  of  flame  spreading  and 
surface  gasification. 

A similar  study  Is  shown  In  Figure  16  for  a larger  crack  (5  = 0.05 
cm).  Because  of  the  large  gap  width,  a very  low  pressure  gradient 
prevails  along  the  crack.  This  causes  the  crack  tip  pressure  to  closely 
track  the  chamber  pressure  before  appreciable  surface  gasification.  For 
this  particular  case,  tpp  Is  more  adequately  defined  by  the  tangent- 
intersection  method.  Similar  to  the  small  crack  results,  the  theoreti- 
cal predictions  for  this  large  crack  reveal  that  a larger  delay  In  t^ 

Is  experienced  when  JPc/dt  Is  reduced. 

The  combined  effects  of  gap  width  and  chamber  pressurization  rate 
are  shown  In  Figure  17.  This  cross-plot  Indicates  that  the  time  required 
to  reach  rapid  pressurization  a-t  the  crack  tip  becomes  shorter  with  a 
decrease  In  the  gap  width  and  with  an  Increase  In  the  rocket  chamber 
pressurization  rate.  For  very  large  cracks,  the  definition  of  t^ 
becomes  less  clear  and  difficult  to  Interpret.  For  extremely  small 
cracks,  there  Is  a critical  gap  width  below  which  a flame  cannot  easily 
propagate;  therefore,  t^  will  Increase  significantly  with  a 
further  reduction  In  gap  width. 

Figures  18  and  19  describe  the  Influence  of  the  rocket  chamber 
pressurization  rate  on  the  flame  front  propagation  Inside  a crack  with  a 
gap  width  of  0.02  and  0.05  cm,  respectively.  It  Is  evident  that  there 
is  a characteristic  difference  between  these  two  sets  of  plots.  For  the 
smaller  crack  (see  Figure  18),  the  strong  convective  flow  enhances  an 
early  attainment  of  Ignition  near  the  crack  entrance.  The  resultant 
propellant  gasification  reinforces  the  convective  flow  field  and  causes 
the  flame  front  to  rapidly  accelerate  along  the  crack.  For  distances 
along  the  crack  less  than  5 cm,  the  rocket  chamber  pressurization  rate 
significantly  affects  the  rate  of  flame  propagation.  However,  for  x > 5 
cm,  the  flame  front  propagation  velocity  (approximately  800-900  m/s)  Is 
dominated  by  Interior  propellant  gasification,  and  the  Influence  of  the 
rocket  chamber  pressurization  Is  diminished;  therefore,  the  slopes  of 


FIGURE  18.  Ini 
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these  curves  are  approximately  the  same.  For  the  larger  crack  (see 
Figure  19) , the  rocket  chamber  pressurization  rate  and  the  local  gasi- 
fication processes  are  both  important  factors  controlling  the  rate  of 
flame  propagation  and,  as  discussed  earlier,  the  pressure  at  the  crack 
tip. 

The  Influence  of  propellant  ignition  temperature  on  the  time  to 
reach  rapid  pressurization  at  the  crack  tip,  tjjp,  is  shown  in  Figure  20 
for  various  rates  of  rocket  chamber  pressurization.  It  is  found  that 
tRp  for  this  small  crack  (6  = 0.02  cm)  diminishes  with  a reduction  in 
the  ignition  temperature.  Also,  the  slope  of  tgp  versus  T^gjj  decreases 
as  the  chamber  pressurization  rate  Increases.  This  implies  that  the 
flame  spreading  and  pressurization  processes  inside  a crack  of  a given 
initial  geometry  become  less  sensitive  to  the  propellant  ignition  temper 
ature  as  the  forced  convection  becomes  stronger. 
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FIGURE  20.  Influence  of  Propellant  Ignition 
Temperature  and  Rocket  Chamber  Pressurization 
Rate  on  the  Time  to  Reach  Rapid  Pressurization 
at  Crack  Tip. 
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SUMMARY  OF  PROGRESS  AND  CONCLUSIONS 


The  experimental  and  theoretical  progress  achieved  at  Pennsylvania 
State  University  In  understanding  the  complicated  combustion  processes 
that  occur  In  a single.  Isolated,  propellant  crack  Is  summarized  In  the 
following. 

1.  A one-dlmenslonal  theoretical  model.  Incorporating  the  Noble- 
Abel  dense  gas  law,  has  been  developed  to  describe  the  transient  combus- 
tion phenomena  In  a propellant  crack.  The  theoretical  model  can  be  used 
to  predict  wave  phemomena,  heat  transfer  from  the  gas  to  the  propellant 
surface  and  associated  thermal  penetration,  flame  propagation  and  resul- 
tant pressurization  at  various  locations  along  the  propellant  cavity. 

2.  Several  Important  results  of  the  theoretical  study  are  listed 
below. 


a.  For  cracks  with  a gap  width  above  the  critical  value 
required  for  gas  penetration,  the  Internal  pressurization 
rate,  pressure  gradient,  and  flame  velocity  in  the  crack 
Increase  as:  the  gap  width  decreases;  the  rocket  chamber 
pressurization  rate  Increases;  and  the  propellant  igni- 
tion temperature  decreases. 

b.  During  the  early  phase  of  flame  spreading,  pressure  wave 
phenomena  exist  and  the  crack  pressurization  rate  is 
controlled  by  the  rocket  chamber  conditions.  At  a later 
time,  the  cavity  pressurization  Is  dominated  by  propel- 
lant combustion. 

c.  The  calculated  flame  front  deceleration  near  the  crack 
tip  agrees  with  experimental  observations  made  separately 
by  Taylor^  and  Belyaev. 6 

3.  A laboratory-size  solid  propellant  combustion  chamber  has  been 
designed,  fabricated,  and  tested.  A series  of  firings  with  single-pore 
propellants  has  been  conducted  for  Initial  system  checkout. 

4.  Preliminary  experimental  data,  obtained  In  wide  cracks  (6  0.1 

cm)  for  relatively  low  chamber  pressurization  rates,  are  In  agreement 
with  theoretical  predictions. 
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NOMENCLATURE 


l! 


2 

Ap  Cross-sectional  area  of  the  crack,  cm 

Apx  Spatial  change  of  cross-sectional  area  of  crack  with  respect  to  the 
axial  distance,  cm 

a Pre-exponential  factor  in  the  nonerosive  burning  rate  law,  ap  , 
(cra/s) (cm2/gf)n 
Bx  Body  force,  gf/g 

b Co- volume,  cm^/g 

c Speed  of  sound,  cm/s;  when  with  subscript,  specific  heat,  cal/g-K 

Cj  Friction  coefficient,  2gTTv,/pu2 

c Specific  heat  at  constant  pressure,  cal/g-K 

d^  Hydraulic  diameter  of  the  crack,  cm 

E Total  stored  energy  (internal  and  kinetic),  cal/g 

g Acceleration  of  gravity,  conversion  factor,  g-cm/gs-s^ 

H Depth  of  the  gap  in  flow  visualization  tests,  cm 

— 2 
h(.  Local  convective  heat-transfer  coefficient,  cal/cm  -s-K 

hj,p  Local  convective  heat-transfer  coefficient  over  the  propellant  sur- 

_ face,  cal/cm^-s-K 

h Local  convective  heat-transfer  coefficient  over  non-propellant  port 

wall,  cal/cm'^-s-K 

hf  Enthalpy  of  combustion  gas  at  adiabatic  flame  temperature,  cal/g 

J Mechanical  equivalent  of  heat,  gf-cm/cal 

Kg  Erosive  burning  constant,  cm^-K/cal 

L Length  of  the  crack,  cm 

M(,  Molecular  weight,  gf/g-mole 

n Pressure  exponent  in  the  non-erosive  burning  rate  law 
Pr  Prandtl  number 

Burning  perimeter,  cm 
^ Wetted  perimeter  of  the  port,  cm 
P Static  pressure,  gf/cm^ 

(or  p) 

R Specific  gas  constant  for  the  combustion  gases,  g^-cm/g-K 
Re  Reynolds  number 

rj^  Burning  rate  of  the  solid  propellant,  including  the  erosive  burning 
contribution,  cm/s 

T Temperature  (without  subscript,  static  gas  temperature),  K 

Taf  Average  film  gas  temperature  (T  + Tps)/2,  K 

Tf  Adiabatic  flame  temperature  of  the  solid  propellant,  K 

Tpi  Initial  propellant  temperature,  K 

Tpg  Propellant  surface  temperature,  K 

'^ws  Nonpropellant  wall  surface  temperature,  K 


t 

u 


x 


XL 

y 


Symb^ 

] 


6 

a 

6 

Y 


e I 

^s  I 

X 

w ■ 


p 


Subs^ 


c 

cri 

i 

ign 

pr 

RP 

tip 


5 


;h  respect  to  the 
ite  law,  ap  , 


I heat,  cal/g-K 


L 


jfJ-s-K 

! propellant  sur- 

l-propellant  port 
\ 

erature,  cal/g 


Law 


cm/g-K 

erosive  burning 
:ure ) , K 
t,  K 


t Time , s 

u Gas  velocity,  cm/s 

Vgf  Velocity  of  propellant  gas  at  the  burning  surface,  cm/s 
Vjj  Coefficient  in  viscosity-temperature  relation 

X Axial  distance  from  propellant  crack  opening,  cm 
XL  Position  at  the  end  of  crack,  cm 

Xp  Axial  distance  along  the  crack  at  which  propellant  begins,  cm 
y Perpendicular  distance  from  the  propellant  surface  into  the  solid,  cm 

Symbols 

6 Gap  width  of  the  crack,  cm 

a Thermal  diffusivity,  cm^/s 

6 Erosive  burning  exponent 

Y Ratio  of  specific  heats 

e Small  number  of  K 

Eg  Surface  roughness,  cm 

X Thermal  conductivity,  cal/cm-s-K 

y Gas  viscosity,  g/cm-s  (poise)  2 

p Density  (without  subscript,  gas  density),  g/cm 

Shear  stress  on  the  port  wall,  g^/cm^ 

Txx  Normal  viscous  stress,  gf/cm^ 

Angle  measured,  in  a counterclockwise  direction,  at  the  lower  side 
of  the  propellant,  degree 


Subscripts 


c Rocket  chamber 

cri  Critical  condition  for  surface  ablation 

i Initial  value 

ign  Ignition  condition 

pr  Propellant 

RP  Rapid  pressurization 

tip  Crack  tip 


f 


L 


mm  m 

' /r,l4fe</  J| 
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